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We examine the dynamical implications of an interaction between some of 
£SJ ' the fluid components of the universe. We consider the combination of three 

matter components, one of which is a perfect fluid and the other two are 
interacting. The interaction term generalizes the cases found in scalar field 
cosmologies with an exponential potential. We find that attracting scaling 
\ solutions are obtained in several regions of parameter space, that oscillating 

behaviour is possible, and that new curvature scaling solutions exist. We 
also discuss the inflationary behaviour of the solutions and present some of 
the constraints on the strength of the coupling, namely those arising from 
O ', nucleosynthesis. 
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I. INTRODUCTION 



There has been a growing appreciation of the importance of the asymptotic behaviour 
of cosmological models 0. Indeed, unless there is a case for a cosmic coincidence 0, the 
features of the dynamics should be associated with some stationary regime which should be 
obtained without fine-tuning of the initial conditions. 

A particular case which has attracted a great deal of interest concerns the possibility 
of obtaining cosmological scaling solutions with self-interacting scalar fields. It was shown 
that the exponential potentials yield the remarkable feature that the dynamics of the scalar 
field self-adjusts to that of matter so that the corresponding energy densities become pro- 
portional 0-0. In these solutions one envisages the interplay between the scalar field and 
matter, instead of focusing on the dynamics of models exclusively dominated by a self- 
interacting scalar field as was the case in most of the models of inflation. It was shown that 
these solutions attract all the other phase space trajectories in the case of flat space mod- 
els and hence provide the late-time asymptotic behaviour for the scalar field cosmologies 
under consideration. This gives, on the one hand, a possible answer for why a non-vanishing 
scalar field does not introduce radical changes with respect to the usual Einstein-de Sitter 
rate of expansion of the universe. On other hand, it may contribute an explanation to the 
difference between the actual density of matter and the critical energy density of the spa- 
tially flat isotropic models. Furthermore, the scalar field component would also fulfill the 
convenient role of delaying the time of matter-radiation equality which would help fitting 
the power spectrum of large scale structure [p|-[T0|]. 

So far the emphasis has been placed on the role of given families of potentials. In 
the literature we find mainly two sorts of potentials underlying the scaling behaviour: the 
exponential potentials and a class of power law potentials |2|,||,|ll[] respectively. However it 
is worth noticing that the solutions corresponding to the latter set of potentials were only 
shown to hold in the regime where the perfect fluid component fully dominates the expansion 
and the energy density of the scalar field is negligible. Regarding the solutions with an 
exponential potential the self-adjustment of the behaviour of the energy density of one or 
more scalar fields with that of matter has been investigated in Friedmann-Robertson- Walker 
(FRW) models [0-0JT^,[13| both with and without curvature, in spatially homogeneous, but 
anisotropic models jnj, and in FRW models in scalar-tensor gravity theories (also referred 



to as non- minimal coupling) fl5[-18 . It was also shown by two of us that every positive 
and monotonous potential which is asymptotically exponential yields a scaling solution as a 
global attractor. 

Recently Billyard and Coley pO have included an interaction term which transfers energy 
from the scalar field to the matter fields. In this work we analyse a related question, namely 
the role of the interaction between two of the components of the universe in promoting the 
scaling behaviour. We show that in the case of the self-interacting scalar field the self-tuning 
is a direct result of the energy transfer between the purely kinetic part and the vacuum-like 
part of the field accomplished by the gradient of the potential. We generalize this case 
by allowing for an interaction between two otherwise barotropic fluids and we examine its 
implications for the cosmological behaviour of homogeneous and isotropic universes. We 
consider a phenomenological interaction term between two of the matter components with 
the form oc if A pfpf , where a, (3 and A are constants which on dimensional grounds have 
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to satisfy X + 2(a + (3 — 1) = 0, and which, as we will show below, extends the case 
found in the scalar field with an exponential potential. Such a kind of energy exchange 
between two components guarantees the global conservation of energy-momentum, that 
is, satisfies the contracted Bianchi identities and is akin to the type of dissipative terms 
considered in single component models where either a bulk viscosity term or matter creation 
are considered [ 23]— 125 1 . We show that the consideration of such an interaction term yields a 
variety of situations where scaling solutions emerge. We produce a full classification of those 
cases where non-trivial scaling solutions (NTSS) are attractors. They include as particular 
case the scalar field cosmologies already referred to. Moreover we show that oscillating 
behaviour is also possible. In this latter case although the mean expansion rate is of a 
power-law type the various matter components oscillate. 

This paper is organized as follows. In the next section we review the scaling behaviour 
associated with scalar field cosmologies. We show how we may view this effect as the result of 
an interaction between the two limit situations where the energy of the scalar field lies in its 
kinetic part and the alternative case where it lies in its vacuum-like part associated with its 
potential. We show how the transfer of energy between these is promoted by the gradient 
of the potential. We also consider scalar field cosmological models within the extended 
framework of non- minimally coupled gravity theories. The Brans-Dicke (BD) scalar-tensor 
theory is remarkable in that in the conformally transformed Einstein frame the original 
coupling between the BD scalar field and the curvature scalar of the space-time is traded into 
a re-defined scalar with an exponential potential which is now coupled to matter. In scalar- 
tensor theories there is then an intrinsic interaction between the conformally transformed 
scalar field energy density and the matter fields energy densities which lies at the root of the 
scaling behaviour of some solutions. In fact we consider the matter fields to be a combination 
of a perfect fluid and a radiation fluid and find a new type of scaling behaviour both in the 
Einstein and Jordan frames. 

In section 3 we generalize the scalar field models to a combination of three matter com- 
ponents where one of them is a perfect fluid and the other two are interacting. We perform 
a qualitative analysis which reveals that non-trivial scaling solutions arise in a number of 
situations and we classify them. 

Finally in section 4 we summarize and discuss in the context of our general model a 
number of issues such as inflationary behaviour, curvature scaling solutions, bounds from 
nucleosynthesis and a phenomenological approach to the decay of massive particles out of 
equilibrium. 



II. SCALAR FIELD COSMOLOGIES 

A. Minimal-coupling theories 

We consider a homogeneous and isotropic flat Friedmann-Robertson- Walker (FRW) uni- 
verses filled with a perfect fluid, characterized by p = (7 — l)p, where < 7 < 2 is a 
dimensionless constant, and a self-interacting scalar field within the framework of Einstein's 
theory of general relativity. The field equations then read 

H 2 = ^ + V^)+p, (1) 
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a . dV(<p) 
a d(p 

p = -3H 1P , (3) 

where the overdots stand for the derivatives with respect to the time t, and H = a/a (we 
use units in which c = 1 = 8n/mp, where m p = (-\/G~)~ l is the Planck mass and G the 
gravitational constant). Another equation which is useful, albeit not independent of the 
former, is 



H 



\ [f + 7P) • (4) 



Defining p K = (p 2 /2 and py = V(<p), the scalar field equation (|2|) may be cast as a system 
of two equations 



p K = -6Hp K - V2 (yj pf 



Pv (5) 



pv = +V2 Opj pfp v , (6) 

which shows how the gradient of the potential promotes an interaction between two limit 
"perfect fluids", namely one associated with the scalar field's kinetic energy, px, that may 
be characterized by 7^ = 2 and hence is akin to a stiff fluid, and the other associated with 
the potential energy, pv, which has a vacuum like character and 71/ = 0. In the absence of 
this interaction, when V is constant (the case of an effective cosmological constant), there is 
a smooth evolution from the early time domination of the perfect fluid with a higher value 



of 7 to the late time domination of the perfect fluid with the lower 7 [^,^]. In this case, 
there are no scaling solutions as the asymptotic solutions occur with the vanishing of one of 
the perfect fluid components. 

Let us further define p c and pd as 

Pc = Pk + Pv = y + V &) ( 7 ) 

Pd = Pk - Pv = y ~ V &) ■ (8) 

which are, respectively, the energy density and the pressure of the self-interacting scalar 
field [26[] in the co-moving observer frame. It then follows from Eqs. (||||) 



Pc = -3H ( Pc + p d ) (9) 
Pd = -3H ( Pc + p d ) - 2V'(<p) Vp7+P~d ■ (10) 



The second equation (|10|) is the evolution equation for the scalar field pressure and we see 
that it is this one which explicitly involves the interaction between pi and p2 as discussed 
above. 

When we take V(ip) to be an exponential potential, i.e., V((p) = V exp (z/<p), where v is 
a constant, and introduce r = In a 3 as the new time variable, equations @ and (0) become 



4 



p' c = -(p c + p d ) (11) 

Pd = '(Pc + Pd)~ {pc - Pd) VPc + Pd ■ (12) 

We further introduce the new variables x and y defined as 

V = ^ (14) 

which correspond to the density parameters associated with p c and with p^, respectively. 
Notice that in previous works on scaling solutions with exponential potentials [|7| |i2l , |l3f it 
became popular to use expansion normalized variables which are square roots of density 
parameters. The only reason for that choice is to avoid the ^fp~K which appears in the 
Eqs. (|5]||). Since, in the next sections, we shall be interested in models where the interaction 
term will involve other powers of the individual energy densities, apart from the power 1/2, 
it is preferable to use the density parameters themselves given their immediate connection 
to observations. 

With the definitions flTBD, (IT3) we obtain 



x' = (7 — l)x — y + xy — (7 — l)x 2 (15) 
y' = -x + (7 - l)y + y 2 - (7 - l)xy + 5 y/x + y (x - y) , (16) 

where 5 = —u/\f?> and we have made use of the fact that the Friedmann constraint equation 
for the flat models, Eq. (|l|), yields Q = 1 — x — y. 

It is straightforward to see that we just need to consider the triangle in the x, y plane 
bounded by the invariant lines x = y, x = —y and x = 1. The vertices of this triangle are 
trivial critical points. The origin (0, 0) represents the vanishing of the scalar field energy 
density and hence of its pressure as well (p c = pd = px = Pv — 0, the universe scales as 
a oc t 2 /( 37 ) and p oc a -37 ). The point (1, 1) represents the case where the kinetic energy of 
the scalar field dominates and thus corresponds to the massless field case (p = pv = 0, the 
universe scales as a oc t 1 ^ 3 and px oc a~ 6 ). Finally, the point (1, —1) represents a spurious 
solution introduced by the change of variables we performed, for which p = px = and 
pv = 3H 2 . Clearly, this is not a solution of the original equations (|l|JlJ|). Notice though 
that in cases where the potential exhibits an underlying non-vanishing vacuum energy this 
critical point corresponds to the late-time domination of the vacuum energy and we have 
the well-known de Sitter exponential behaviour |2T||2"7|] . 



Defining e = 7 — 1 it is easy to verify that there are critical points that correspond to 
x = 1 or y = ex. The first one, located at (x = 1, y = 5 2 — 1), is a scalar field dominated 
solution (p = 0) and only exists for 5 2 /2 < 1 which translates into v 2 < 6. This solution 
is an attractor and it corresponds to the well-known power-law inflationary solutions when 
< 2 [29|-31,14 . The second one is given by (x = (1 + e)/5 2 ,y = e(l + e)/5 2 ) is the 



non-trivial scaling solutions found by ||,||] and also studied by and by 0. It only exists 
for v 2 > 37 since x < 1 QJ^] . Linear stability analysis permits to distinguish two topological 
behaviours associated with this fixed point. It is a stable node if v 2 < 24^/(9^ — 2) and a 
stable focus otherwise. 
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B. Non-minimally coupled scalar fields 



Scaling solutions in non-minimally coupled theories have also been found in the literature 



15 — 18 1 . These theories, which can be formulated as general scalar-tensor gravity theories, 



are based on the Lagrangian p2| 

U = $R - g ab <5>,a<5>' b + 2*7($) + 16vrL m , (17) 

where R is the usual Ricci curvature scalar of a spacetime endowed with the metric g a b, $ is a 
scalar field, u;($) is a dimensionless coupling function, U(Q>) is a function of $, and C m is the 
Lagrangian for the matter fields. They provide a most natural generalization of Einstein's 
general relativity (GR), and their investigation enables a generic, model- independent ap- 
proach to the main features and cosmological implications of the unification schemes which 
involve extra-dimensions. 

A distinctive feature of these theories is the coupling of the dynamical scalar field $ 
to the scalar curvature R, which implies that the gravitational constant is now a function 
of $, in fact G = The archetypal of these theories is Brans-Dicke theory (BD), for 

which is constant The Lagrangian (|TT| ) corresponds to the so-called Jordan frame, 
in which the matter fields satisfy the equivalence principle (hence their energy-momentum 
tensor satisfies VbT afc = 0). This means that the L m terms do not explicitly involve the 
scalar field By means of an appropriate conformal transformation of the space-time 
metric g a b to the so-called Einstein frame, 

gab -> 9ab = (^/^*)gab , (18) 

where is a constant allowing the normalization of Newton's constant in the latter frame, 
we recover a minimally coupled theory where the coupling of $ to the curvature is traded 



into a coupling of the redefined scalar field with the matter fields. In fact $ — > ip p7|j38 
such that 



dln$ 16vr , , 



where a = (^2u(<p) + 3) 1 . The field equations for the flat FRW universes with a perfect 
fluid reduce then to the simple form 



•2 

a 2 



. 2 

^ + d~ 3 m^) + v{<p) 



(20) 



^ + ^ = -a-3« + ^M, (21) 
a dtp d^p 

P = -3— P + a~^— ^0 22 
a dip 

where the overdots stand for the derivatives with respect to the conformally transformed 
time t, V(tp) = $1 £/($(y?))/(87r$ 2 (^)), M(ip) = pi ($(v?)/$*) _(2_37/2) and fi is the constant 
defined by pi = pa 37 which fixes the initial conditions for the energy-density of the perfect 
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fluid (note that for the sake of making clear the effect of the conformal transformation on 
the coupling we have written the latter equations with the 87r/$* in spite of of our choice of 
units; in what follows we shall again set it equal to 1). From Eq. ([Zl]) and the definition of 
M(ip) it is apparent that if the matter sources are radiation fields (the case in which M = p 
is constant) or, alternatively, vacuum (the case in which p = and hence M — 0) the 
coupling vanishes, which translates the fact that the latter cases are conformally invariant. 

It is a simple matter to see from Eq. (|TJ|) that in the BD case, where uj (and hence a) 
is constant, $ oc exp (^y/2aipj so that M is exponential and so is V if the original potential 
[/($) is a power-law in $. This means that the coupling encapsulated in the M(tp) function 
amounts in this case to a modification of the constant coefficients of the dynamical system 
associated with the general relativistic case yielding scaling solutions for suitable values of 



the parameters ||i5|--|r8|j . 



Here we shall restrict ourselves to the original BD theory where, besides having a constant 
u (or equivalently a), we also have a vanishing U(Q) in Eq. (|H]). We will continue to focus 
on the flat Friedmann-Robertson- Walker model, but instead of just having a perfect fluid 
we will also consider the simultaneous presence of radiation (hereafter we shall drop the 
tildes referring to the conformal frame quantities). Denoting the perfect fluid (respectively, 
radiation) energy density and pressure by p m andp m (respectively, p r &ndp r ), and the energy 
density and the pressure of the redefined scalar field ip by p^ and p^, the field equations are 

3H 2 = p r + P v + Pm , (23) 

dH 1 

~fa = ~2 + ^) + (P<P + Pv) + (P m + P™)) ' ( 24 ) 

and the evolution equations for the energy densities of the three fluids are given by 

^ = -±H Pr (25) 

IT = - 3H{pm + Pm) + Ti^+I {3pm - pm) lt> (26) 
^ = ~ mp * - {3pm - prn) \t ' (27) 

rendering explicit the interaction that now leads to a transfer of energy between the matter 
perfect fluid and the stiff fluid (0 2 /2) associated with the massless scalar field tp. Since 
Pm = (7 — l)Pm an d given the definitions of p m and of p v , we see that the interaction term 
has the same form as in the case of the minimally coupled scalar field with an exponential 
potential. Indeed, we may write for each component an equation of state of the form 
Pm = (7 — l)p m — Pint and Pip = Pip + Pint, where the interaction term can be cast as 

Pmt = (3H)- 1 J— (4 - 3 7 )p mv //^ • (28) 

There are two major differences with regard to the minimally coupled case. On one hand, 
the dimensionless constant 5 that measures the strength of the interaction depends on 7 
according with 

1 4-37 



2-7 



a . (29) 
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On the other hand, the 7-index of the non-interacting fluid is now fixed to take the value 
4/3 of radiation and 71 remains equal to 2, and 72 now takes the value of the free parameter 
7. Therefore taking into account these differences the dynamical system ( ]T5|JT6| ) adopts the 
form 

, (2 - 3 7 ) , (2 - 3 7 ) 2 

x =— -x — y + xy ; -x 30 

3(2-7) y y 3(2-7) v ; 

(2 - 3 7 ) , 2 (2 - 3 7 ) 



V = -x + - y + y - _ ^ xy + 5 ^x + y (x - y) . (31) 

It is a simple matter to see that as before we obtain two NTSS. One at the point (x = l,y = 
5 2 /2 — l) which corresponds to a late-time domination of the coupled components, the scalar 
field and the 7-fluid. Radiation is depleted and the universe scales as power-law a oc t 2 /( 3r ) 
where T = 7 + S 2 /2 (1 — 7/2) is defined by the position of the fixed point. Another NTSS 
exists at the point (x = 2/(v / 35«), y = (2 — 37)/(3(2 — 7)) x where the three components are 
simultaneously present and the energy densities of the two that are linked by the interaction 
have adjusted to the radiation behaviour. Thus the scale factor evolves as a oc 

t l/2 

as m 

a universe dominated by radiation, but with the major difference that both the BD-scalar 
fluid and the 7-fluid remain in relevant proportions (this scaling solution corresponds in fact 
to the Crm fixed point in the paper of Amendola f4TJ). 

If we reverse the conformal transformation we are able to recover the exact solutions in 
the Jordan frame (which is commonly taken to be the physical frame). The case of the fixed 



point on the line x = 1 corresponds to the power-law solutions found by Nariai [41] and 



includes the dust solution (7 = 1) derived by Brans and Dicke [36]. The other non-trivial 



scaling solution corresponds again to a power-law in the Jordan frame, but differs from the 
single-fluid exact solutions found in the literature f42| , |4"3"l| . For this we have 



a(t) oc t ±V ~ 2 ' 2 (32) 
$(t) oc t ±2 ^-^ , (33) 

where a and t are the scale factor and the time coordinate in the Jordan frame. This exact 
power-law solution is new and it differs from the corresponding behaviour of the radiation 
BD models for which a oc 



III. COSMOLOGIES WITH INTERACTING FLUIDS 

In this section we shall construct a general model which includes the two previous classes 
of models as particular cases. We assume the matter content of the universe to be the 
combination of three components. One that evolves without interaction with the other two 
and which satisfies the usual barotropic p = (7 — l)p law (in the particular case of the 
non-minimal coupling models of the last sub-section it would correspond to the radiation 
component). The two remaining components are mutually coupled by an interaction term 
of the type 

p i = (>y i -l)p i ±T,H*fip$ , (34) 
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where i,j = 1,2, and a, (3, A are constants, and r] is a dimensionless constant. The ± sign 
of the interaction term means that if it takes one of the signs for one of the components it 
necessarily takes the opposite sign for the other, thus ensuring the overall conservation of the 
energy-momentum tensor V&(T^+T^) = as required by the Bianchi contracted identities. 
Moreover, using dimensional considerations, it is a simple matter to verify that a, (3 and A 
have to satisfy A + 2 (a + (3 — 1) = for the interaction to be a pressure. The dependence 
of the interaction term on the products of powers of the densities of the particle species 
reflects, on one hand, the fact that one expects that it should be proportional to the number 
of collisions between the particles of each species and thus on the product of their number 
densities, and, on the other hand, the fact that in the collisions the energies of the particles 
will be shifted. Furthermore the factor depending on H accounts for the characteristic time 
of the individual interactions. Since a detailed relativistic kinetic model of the interaction 
between two fluids is missing, we keep our model as general as possible by using the two 
free parameters of the set a, (3 and A to allow for unknown aspects of the interaction. The 
model under consideration aims at a phenomenological description of transient phases of the 
universe which arise when some of the material components are not in thermal equilibrium. 
One must be wary that this is in fact a likely situation as the expansion of the universe is 
permanently trying to pull the matter fields out of equilibrium |2"I,S|. In this sense the 



type of interaction generalizes those found in the literature when viscous pressure and matter 
creation terms are considered in the literature on dissipative isotropic models [22- |2"4 , |21 ,25 . 
Since our main purpose is to produce a classification of the cases which lead to scaling 
solutions, where two or more components self-adjust their behaviour so that their energy- 
densities scale with the same rate, we shall be essentially concerned with the dynamical 
aspects of this model, rather than with the specifics of particular models. Nevertheless we 
emphasize that the post-inflationary reheating period, the decay of massive particle species 
into lighter ones and the general situations of decoupling of particle species provide examples 
of the situations where the present model may be applied. We will address some of these 
examples in the section IV. 

We shall denote the energy density of the first component p x and assume that it has a 
barotropic equation of state characterized by a constant adiabatic index •y x which we shall 
leave free parameter. 

The field equations read 

3H 2 = P x+Pi+p2 (35) 

H = ~\ [{Px + Px) + (Pi + Pi) + (p 2 + P2)] (36) 
p x = -3H(p x +p x ) (37) 
y £p i = -3H y £(jH+Pi). (38) 

i i 

Here pj, pj, with j — 1,2 are, respectively, the energy density and the pressure of the j th 
matter component measured by a comoving observer. 
We define 

Pc = Pi + p2 (39) 
Pd = Pi - P2 (40) 
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7i + 72 ,,,, 
7c = — g ^ 

7i - 72 ,„ , 
7d = — ~ 2 — (42) 

so that the field equations become 

3H 2 = p x + Pc (43) 
H = - - [7xPx + 7c Pc + 7d Pd] (44) 

px = -3H-f X px (45) 
p c = -3# (7 C p c + 7 d p d ) (46) 
p d = -3# ( ldPc + lcPd ) + 37/F A+1 2 1 - a -^(p c + p d ) a (p c - p d )^ . (47) 

It is also convenient to introduce the density parameters 

* - im < 48 > 

» = ^ (49) 

as dynamical variables and the new time variable r = In (a/a G ) 3 . Then, from Eq. (pE3"D, 
px/3H 2 = 1 — x and the field equations reduce to the following dynamical system 

x' = - ((7c - lx)x + Jay) {x - 1) (50) 
2/' = -7d* - 7c2/ + V (7x(l - a:) + 7c* + l d y) + Q) rj H x+2a+2 ^~ 2 (x + y) a (x - yf , (51) 

where the prime denotes differentiation with respect to the r variable. 

It is remarkable that the dimensional relation that must be satisfied by the parameters, 
X + 2(a + /3 — 1) = 0, is precisely the one that renders the above system autonomous. From 
(|43|), (|48|) and (f49|), we have the restrictions < x < 1, \y\ < x which define the phase space 
of the system, a triangle in the x, y plane bounded by the invariant lines x = y, x = —y 
and x = 1. Using the symmetry (y,jd, ®,(3) — > (—y, — 7d, A oc) we may restrict ourselves to 
the when 7^ > 0, i.e., 71 > 72. Rescaling time through the factor 7^ and defining the new 

parameters e = — — — , 5 = r L^— . equations (pDI) and (pll) may be simplified and 

Id Id 
rewritten in a more compact form as 

x — (1 — x)(—y + ex) 

y = —x + ey + y(y — ex) + <5(x + y) a (x — y) 13 . (52) 

The parameter 5 measures the strength of the interaction, and |e| < 1 indicates that 7x 

lies within the interval (71,72). The cases of the previous section correspond to a = 1/2, 

(3 = 1. In the minimally coupled scalar field cosmologies, e = 7 — 1, whereas in the BD case 
, _ 2-37 

6 ~ 3(7-2)- 

The points (0, 0), (1, 1) and (1, —1) are always equilibrium points for this system. They 
correspond to the trivial scaling solutions which exist when rj = 0. We shall say that an 
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equilibrium point of (|5^) is a non trivial scaling solution (NTSS) if it is none of the above 
and furthermore it is stable. Our study of system (|52|) will be directed towards the search 
for non trivial scaling solutions. The additional symmetry (y, t, e, a, (3) — > (—y, —t, —e, (3, a) 
is useful to reduce the number of different possible cases. The techniques we employ are the 
standard methods of qualitative theory of planar systems, assisted, in the degenerate cases, 
by numerical integration. 



A. Dynamics on x = 1 



When px = 0, the interplay between the two interacting fluids is described by the single 
equation 

y = (yZ-l) + 5(l + yr(l-yf. (53) 

It is straightforward to check that NTSS on x = 1 show up in the three cases depicted in 
Figure 1. Denoting generically by (x,y) the coordinates of any of these fixed points, we see 
from Eqs. (|43|), ([46|), (fl8D and (^) that the scaling solution is characterized by a oc t 2 /( 3r ) 
where T = j c + j d y. 



B. Global dynamics for |e| > 1 



When |e| > 1 there are no equilibria of (52) outside the boundary of the phase space. 
For e < —1, i.e., 7x < 72 < 7i, the origin is a global attractor. For e > 1, the line 
x = 1 becomes the global attractor, so that the stable equilibria described in the preceding 
subsection become NTSS for the full system, which behaves as shown in Figure 2. 



When the non-interacting matter component behaves as one of the two interacting fluids, 
we have the degenerate cases e = 1 (jx — 7i) an d e = — 1 (jx = 72)- These cases also 
provide NTSS on the lines y = x and y = —x, respectively, as shown in Figure 3 for the 
cases when other equilibria on x = 1 coexist. In the cases of Figures 3. a) and 3.b), that is, 
when 7x = 72 and a > 1, px would have been negligible in the past, while in the case of 
Figure 3.c), that is when ^ x = 7i and /3 < 1, px would have dominated in the past. 



C. The case |e| < 1 and a,/3 > 1 



When (1 + e) Q " 1 (l — e)^ 1 < 1/5, there are no equilibria in the interior of the phase 
space, and the dynamics is either trivial or as shown in Figure 4. a) according to whether 
{ a+ 2 p-2 ) a+l3 ~ 2 ( a ~ l) a ~ 1 (/5 — 1)^ _1 is smaller or greater than 1/5, respectively. When (1 + 
t) a ' l (l - ef- 1 > 1/5 we have always (^Z2) a+t3 ~ 2 (a - l) 0-1 ^ - l)" -1 > 1/5, since the 
first member of the preceding inequality is precisely the maximum of (1 + e) a_1 (l — e)' 3-1 
as a function of e, attained at €m = Q °^ 2 . In this case, there exists another equilibrium in 
the interior of the phase space, whose stability is determined by the sign of eu — e- Thus, 
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we have an additional NTSS in the cases of Figure 4.b) and 4.c). The coordinates of the 
new solution ), where x satisfies 

l/5 = x a+fi -\l + t) a -\l-tf- 1 . (54) 

At this scaling solution and, in fact, at every scaling solution in the interior of the phase 
space domain, the universe scales as a oc t 2 ^ 3/yx \ 

In the case of Figure 4.b), we have the new feature of a neutrally stable NTSS. In the 
case of Figure 4.b), the basin of attraction of the NTSS shares the phase space with the 
basin of attraction of (1,-1). 

D. The case |e| < 1, a,/3 < 1 and a + j3 > 1 

As above, an equilibrium in the interior of the phase space exists only when (l + e) 1 ~ a (l — 
e) 1 ^ < 5. For (l + e) 1 ~ a (l -e) 1 ^ > 5 we have always (^l^f+^a-l)"- 1 ^-!)^ 1 > 5 
since, as before, this is the maximum value of the function (1 + e) 1_a (l — e) 1_/3 , attained 
at tu = (0 — a)/(2 — a — (3). The dynamics in this case is as shown in Figure 5. a). With 
respect to Figure 2.c), the change in e has turned the NTSS into the only attractor. 

When (1 + e) 1 - Q (l - ef~ p < 5, we may have either (1 + e M ) 1 " a (l - e^) 1 ^ < 5 or 
(1 + ejvf) 1_a (l — tM) 1 ^ 13 > <5- In the first case, there are no NTSS on the boundary of the 
phase space, and we have at (x, ex), x given by ( |54"D a neutrally stable (resp. asymptotically 
stable) NTSS when e = €m (resp. e < €m), see Figures 5.b) and 5.c). Notice that, in both 
situations, (0, 0) is an attractor for a set of positive measure of initial conditions. In the 
second case, the dynamics is as shown in Figures 5.d) and 5.e), according to whether e < 6m 
or e > €m- In both situations, the NTSS is a global attractor. 

E. The case |e| < 1, a, (3 < 1 and a + (3 < I 

In this case, we shall have an equilibrium in the interior of the phase space at (x, ex with 
x given by ([54]) whenever (1 + e) 1_Q (l — e) 1_/3 > 5. When this happens, there are always 
equilibria on the line x = 1, which falls into the case of Figure l.c), and the additional 
equilibrium is a saddle. The dynamics is shown in Figure 6. a). With respect to Figure 2.c), 
the NTSS survives as a partial attractor. 

When (1 + e) 1_a (l — e) 1_/3 < 5 there are no equilibria in the interior of the phase space 
and we may still have either (^z^) 2 """^! - a) 1 ~ a (l - /3) 1 " /3 > 5 or (^b^) 2 " "^! - 
ct) 1-a (l — /3) 1_/3 < 5. Only the first case may provide NTSS, when, moreover, e > 0, and 
the dynamics is then as shown in Figure 6.b). With respect to Figure 2.c), the dynamics is 
only slightly changed, and the NTSS. 

F. The case |e| < 1 and a < 1 < f3 or a < 1 < (3 

An equilibrium in the interior of the phase space exists if and only if (l+e) a ~ 1 (l — e) 13 ^ 1 > 
1/5. In this case, it is easy to check that the equilibrium on x = 1 is a saddle, and that 
the additional equilibrium is a sink and hence a NTSS. In the complementary case, there 
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can only be non trivial scaling solutions when the dynamics on x = 1 is as in Figure l.b). 
Then, the NTSS on x = 1 is a sink for the full system, and a global attractor. The dynamics 
is represented on Figure 7. Finally, it is easy to check using the symmetries that the case 
ft < 1 < a provides no NTSS. 



G. The case |e| < 1, a = ft = 1 and 5 > 1 

In this case, we have also a neutrally stable NTSS in the interior of the phase space. The 
dynamics is shown in Figure 8. With respect to Figure 5.b), the behaviour is similar, except 
that (0, 0) is now a saddle so that, in the present case, there are no attractors. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have produced a qualitative study of the dynamics of the flat FRW model 
with three matter components one of which is barotropic and the other two are coupled, 
exchanging energy. The consideration of the flat model is motivated by observational data 
that indicate that the universe is flat p5|-^1 . Since these results also suggest that the density 
parameter of matter (including baryonic matter and cold dark matter) only accounts for 
Q m ~ 0.35, some additional form of dark energy should also intervene. Irrespective of the 
particular proportions taken by each of the major material components the mere fact that 
any other components might be present apart from the usual barotropic radiation and matter 
fluids of the standard cosmological model [f48"l , f49| raises the worry that the usual expansion 



rates during the radiation and Einstein-de Sitter epochs of the standard model would be 
significantly disrupted. Attracting scaling solutions provide a possible way of reconciling 
the presence of several matter components with the standard model as the various energy 
densities at stake all scale with the same rate as that of the barotropic fluid if the latter 
is non-vanishing. For the general type of interaction we considered in this paper, we have 
presented all the classes of parameter values for which system (|52| ) exhibits non trivial scaling 
solutions. Those which correspond to equilibria outside the line x — 1, that is, for which 
px 7^ 0, are especially interesting from the physical point of view. These are the cases of 
Figures 3, 4.b) and c), 5.b), c), d) and e), 7.a) and 8. Among these, the cases of Figures 4.b), 
5.b) and 8 are also remarkable. They correspond to situations where for a set of positive 
(or even full) measure of initial conditions in phase space the solutions are periodic, and so 
the relative abundance of each of the three fluids oscillates. We would also like to remark 
that this class of models opens the possibility of having simultaneously scaling and periods 
of inflationary behaviour. 

A. Inflation 

In fact, it is easy to see that the condition for inflationary behaviour (d > 0) is 

»<« + ^*, (55) 
old 
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and thus defines a half-plane in the x, y plane. This condition does not depend on the values 
of the parameters A, a, (3 and S, it just depends on 7 x, 71 and 72, and adjusting the latter 
the triangle of the physical phase space may have partial or full overlap with the inflationary 
half-plane. When the boundary line y*(x) = ex + (2 — 3 / yx)/(5 / ld) lies in the y > half- 
plane and intersects the x = 1 vertical line above the point (1, 1) all the trajectories are 
inflationary. This situation arises when 7^ < 2/3 and 71 < 2/3 (recall that we are assuming 
72 < 7i)- The opposite situation arises when the boundary line y*{x) lies in the y < 
half-plane and intersects the x = 1 vertical line below the point (1,-1) in which case no 
trajectory is inflationary. This happens for 7^ > 2/3 and 72 > 2/3. In the intermediate 
situations, namely when 72 < 2/3 < 71, the inflationary region corresponds to the portion of 
the phase space triangle below the y*(x) line. The intersections of this line with the triangle 
are at the points 

Old 

1 2-3 7x 1 2-3 7x , , 2 

1 5 , z 5 ) when 7x < - , 

1 - e 3j d 1 - e 3j d 3 

1 3 7x -2 1 2-3 7x , , 2 

■) when 7x > - • 



+ e 3j d 1 + e 37 d ' " " 3 



Let us consider what happens in the scalar field models previously discussed in Section 
II. In the minimal coupling case the inflationary solutions lie below the y*(x) = (7 — l)x + 
(2 — 3 7 )/3 line. This line intersects the x = 1 frontier of the triangle at y — —1/3 regardless 
of its slope defined by the value of 7. For 7 = 1, that is, for dust, the line is horizontal and 
the inflationary region is a triangle. As we consider smaller values of 7, the slope of the line 
y*(x) becomes increasingly negative and in the limit case of 7 = it becomes —1. Notice 
that this limit value of 7 does not yield NTSS solutions in the interior of the phase-space 
domain. It corresponds to a cosmological constant and it is equivalent to having a non- 
vanishing vacuum energy in the exponential potential. The converse happens as we consider 
values of 7 > 1. The slope increases up to a maximum value when the limit 7 = 2 is chosen 
and the boundary line is parallel to y — x. So there is always an inflationary region in the 
phase diagram of the models and this means that portions of the trajectories approaching the 
attractors will exhibit inflationary transients. Let us consider now the question of whether 
the scaling solutions fall within those regions. In the v 2 < 37 < 6 case, the only fixed 
point lies on the x = 1 vertical line. When v 2 < 2 its y coordinate satisfies y < —1/3 and 
hence falls within the inflationary region. In the v 2 > 37 case, the stable scaling solution 
x = 37/z/ 2 , y = 37(7 — l)/^ 2 belongs to the inflationary region only if 7 < 2/3. This 
means that we may have NTSS which exhibit inflationary behaviour. However, as the most 
interesting models from the viewpoint of the late time behaviour of the universe are those 
for which the perfect fluid has 7 > 1, namely 7 = 1,4/3, the remarkable issue is that a 
non-negligible set of solutions naturally undergo a finite period of inflation before reaching 
the attractor. 

The models leading to oscillatory behaviour are also interesting in what concerns infla- 
tion. Consider the model associated with the figure 8, that is, the model characterized by 
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|e| < 1, a = (3 = 1 and 5 > 1 . The NTSS is located at x = 1/5, y = e/5 and it corresponds 
to an expansion that tracks the non-interacting perfect fluid so that a oc 

t 2/(3 7 x)_ Therefore it 

is immediate to see that this NTSS falls within the inflationary region of the phase-diagram 
if 7x < 2/3 and off it otherwise. Taking into consideration what was expounded at the be- 
ginning of this subsection, in order for the inflationary half-plane to overlap the phase-plane 
triangle one requires that 72 be smaller than 2/3. Assuming 7x > 2/3 and 72 < 2/3 we 
have then that the oscillatory trajectories beyond certain radius from the fixed point will 
undergo cyclic periods of inflation since they periodically cross the inflationary region of the 
phase-plane. Notice that since the model under consideration corresponds to a non-linear 
oscillator the periods associated with the trajectories increase from the immediate (and ex- 
tremely small) neighbourhood of the fixed point, where the linear approximation holds true 
and T ~ 27r/e(l — 1/5)JJ\ — e 2 )/(5 — 1), to the regions close to the limits of the triangle 
where the period becomes infinite. 

B. Curvature scaling solutions 

An interesting case which also emerges from our analysis regards the effect of a non- 
vanishing spatial curvature in models with two interacting fluids. Indeed taking the partic- 
ular case where jx = 2/3 for the non-interacting perfect fluid from eq. ( |37|) we see that this 
is equivalent to having a term p oc cT 2 in the Friedmann equation. Conversely, it is easy 
to verify that the usual curvature term taken as pk = —3k/ a 2 satisfies Eqs. (^5]), Q36|) , ( |37|) 
with 7 = 2/3. Thus the case of two coupled fluids in k = — 1 models falls within the scope 
of our study. 

From the definition of e we have then e = (2 — 37J / (Sjd) and so, according to our results 
of Section III, curvature scaling solutions, as defined in ||13|| , exist when 72 < 2/3 < 71. In 
this case the scale factor evolves as a oc t and both p\ and p2 self-adjust to the p\, pi oc a~ 2 
behaviour of the curvature term. If both 71 and 72 are greater than 2/3 the origin is a 
global attractor and we recover the usual asymptotic behaviour found when we have two 
non-interacting perfect fluids. The curvature term eventually dominates and we have the 
vanishing of p\ and of p2- This solution corresponds to well-known Milne universe. 

In the limit case where e = — 1, and hence 72 = 2/3, the curvature term dominates in the 
future, but only a limited set of solutions corresponds to the depletion of both the pi and 
the P2 components. Indeed almost all solutions end up in the x = —y border line and they 
correspond to solutions with a oc t which occur with the depletion of the pi component. In 
this case there are no inflationary solutions in the sense that the scale factor does not evolve 
with a power greater than 1 (we are in the limit coasting model). In the alternative limit 
case, when e = 1 and hence 71 = 2/3, the curvature was the dominating component in the 
past and the cosmological models evolve either towards scaling solutions characterized by 
the depletion of the P2 component or towards scaling solutions where the curvature vanishes. 
In this case all the solutions are inflationary. 
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C. Nucleosynthesis 



It was pointed out in previous works on scaling solutions |3]-|7| that the most stringent 
bounds on the admissible densities of the components that are present in addition to the 



usual perfect fluid are set by the primordial nucleosynthesis of light elements |49]. Since the 
perfect fluid is radiation, the attractor solution is characterized by the usual of expansion 
a oc t 1/2 and so any deviations from the standard model yields of the light elements are a 
consequence of the number of degrees of freedom N(t nuc ) which are due to the extra matter 
components 0. The limits that beset this number can be translated into a permitted range 
of energy density associated to the additional matter contributions which is Q extra <0.13 — 
0.2 g. ' 

In the present case VL extra = x so that we have the following bounds on 5 at the attractor 
scaling solution of Eq. fl54|) 

013^ (56) 

when a + (3 > 1, and 

<5<(l + e) 1 - a (l-e) 1 - /3 0.13 1 - Q - /3 (57) 
when a + (3 < 1. In these expressions e = (4 — /^>ld) since 7x = 4/3. 



D. Decay of massive particles 

We now briefly consider the possibility of applying the present model to a transient regime 
during the early universe when two particle species interact. For instance the decay of some 
massive particle species into a lighter one occurring out of equilibrium. This question has 



been analysed in the literature (for a review see Chapter 5 of [49 1 and references therein) and 
usually it is assumed that the massive particles decay into relativistic particles that rapidly 
thermalize so that, on one hand one may consider them as being a part of the radiation 
component, and on the other hand one does not have to consider reverse processes. 

The model envisaged in this work enables one to relax the assumption of thermal equi- 
librium of the lighter species. Following [49| if we denote by if) the decaying massive particle 



species, the relevant equations which are usually adopted are 

p^ + ZHp^, = -r (58) 
p r + AHp r = r pf (59) 

where p^ is the energy density of the if) species, p r is the energy density of the radiation fluid 
which includes the thermalized daughter products of the decays, and the time at which the 
decays take place is given by r ~ t ~ H -1 . 

If we apply our model it becomes possible to consider the transient stage during which 
the lighter particle species energy density increases due to the decays of the ^>'s and is not 
yet in thermal equilibrium with the radiation. Since there is no well-established thermo- 
kinetic prescription for this situation the purpose of this analysis is mainly illustrative and 
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a more detailed investigation of the specifics of the process is left to a future work. From 
the viewpoint of our model the usual treatment found in the literature may be associated 
with taking p x as being the radiation perfect fluid, to be characterized by 7 2 = 1 and the 
daughter products to be described by pi with 71 = 4/3. Then we take r ~ r/H x+1 Pipf" 1 
in Eq. (|59| ) where the parameters are left free. Thus we have a model characterized by 
e = 1 and we immediately know that the dynamics corresponds to the phase-plane of the 
Fig. 3 c) where only the trajectories above the separatrix that connects the (0, 0) singular 
point to the saddle scaling solution on the x = 1 border of the phase diagram evolve towards 
equilibrium with the radiation component. These solutions correspond to a depletion of the 
Px component and, hence, of the decaying massive particles. 

If we relax the assumption that the daughter products are in equilibrium with radiation 
and we keep the 71 parameter free, we may still have a non-trivial scaling solution in the 
interior of the phase-plane domain, for which both of the components interacting through 
the decays swiftly adjust themselves to track the perfect fluid behaviour of radiation. We 
would have then their thermalization and there is no depletion of the massive particle species 
in this case. For this to happen the lighter species must be characterized by a 71 > 4/3 
in order to satisfy the condition |e| < 1. According with Eq. (|54|) the combined density 
parameter settles at a value determined by 71, since in the present case it completely defines 
e (we have e = (5 — 37i)/(3(7i — 1)), and also by the values taken for a, f3 and 5 (the latter 
is given by 5 = 3(3H /2) a+l3 ~ 2 p^ 1 pi _/3 r _1 ). However as presented in our study additional 
conditions must also be met in order to have an attracting NTSS. The relevant cases are those 
represented in Figs. 4 c), 5 c) and d) and 7 a) and a common condition 5 > (l+e) 1_a (l— e) 1_/3 
which translates into 




(60) 



If one selects particular values of the remaining parameters this means yet another constraint 
on S. 

The main point to emphasize is then that there exists the possibility within the framework 
of a flat model of having a dynamical thermalization of the three components and, hence, 
the energy densities of the particle species adding to the radiation energy density. In the 
conventional models the relics of the ?/>'s do not contribute to the latter and are subject to 
the so-called Lee- Weinberg bound p| . 



E. Conclusion 

To conclude we want to stress that our results reveal how the consideration of nonlinear 
interactions between some of the matter components of the universe allows for the emer- 
gence of a variety of phenomena of which the scaling solutions and oscillatory behaviour 
are remarkable examples. In this sense, our study resembles the approach of Ref. [ 50f 



m 



their analysis of the dynamics of models with bulk viscosity. Apart from producing the 
classification of the models in parameter space in terms of the qualitative behaviour of the 
solutions, we have also addressed some of the cases that can be singled out not only from the 
point of view of their dynamics, but also from the perspectives of possible application to the 
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thermal physics of the universe. From this more physical standing, we would like to stress 
the following aspects. First, our models generalize the scalar field models yielding scaling 
behaviour thus providing a different physical setting in which the interesting properties of 
scaling solutions may be obtained. Second, the classification scheme provides the guide- 
lines to be followed in the search for the possible causes of a given phenomenology. Third, 
the detection of new phenomena or the need to re-evaluate observational limits regarding, 
for instance, relic abundances, nucleosynthesis, dark matter, may be reexamined with the 
help of the classification presented here. However, needless to say that the consideration 
of any specific model requires a detailed kinetic analysis to provide a solid justification to 
the particular interaction model. Finally, our study has lead naturally to new results on 
the existence of curvature scaling solutions (sub-sect. IV-B) and on the conditions for finite 
inflationary periods (sub-sect. IV-A). 
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FIGURES 



a) b) c) 



FIG. 1. Non trivial scaling solutions on the invariant line x = 1. a) a, (3 > 1, 



{^T=2) a+I3 ~ 2 {a - ly- 1 ^ - If' 1 > \. b) a < 1 < (3 or a = 1, (3 > 1 and 2^ > ±, 



or 



a < 1, /? = 1 and 2^ < J. c) a, (3 < 1, (^-^/^(i _ a )a-i (1 _ ^,9-1 < i . 



a) b) c) 

FIG. 2. Global dynamics for e > 1. a) a,/3 > 1, (^|r2) a+/3 ~ 2 (a - l) a ~H(3 - I)' 3-1 > y- 

b) a < 1 < (3 or a = 1, (3 > 1 and 2 /3 ~ 1 > ±, or a < 1, /? = 1 and T-- 1 < \. c) a,(3 < 1, 
(^^(i-aja-^i-^-^i. 



a) b) c) 

FIG. 3. Non trivial scaling solutions for |e| = 1. a) a, (3 as in Figure 2. a) and e = — 1. b) e = — 1 
and (3<l<aoi/3 = l, a > 1 and 2 a_1 > |, or (3 < 1, a = 1 and 2 /3 ~ 1 < j. c) a, /3 as in Figure 
2.c) and e = 1. 
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FIG. 4. Non trivial scaling solutions for |e| < 1 and a,/3 > 1. a) (1 + e) a 1 (1 - e) 13 ^ 1 < 1/5 
and (^W^) a+ ^ 2 ( a ~ l) Q_1 (/3 ~ > 1 / S - b ) C 1 + e) Q_1 (l " e)' 3 " 1 > 1/* and e = e M - c) 

(1 + e)""^! - ef- 1 > 1/5 and e > e M - 




a) b) c) 



d) e) 

FIG. 5. Non trivial scaling solutions for |e| < 1, a, (5 < 1, and a + /3 > 1 
. a)(l + ey- a (l - e) 1 "" > 5. b),c) (1 + e) 1 - Q (l - e) 1 ^ < 5 and ( 1 -|-^) 2 -«-/ 3 (l - 
- /?) 1 ~ /3 < 5. b) e = e M - c) e < e M . d), e) (1 + e) 1 " a (l - e) 1 "" < 5 and 
(db) 2_a_/3 (l - - > 5- d) e < e M . e) e > e M - 



a) b) 

FIG. 6. Non trivial scaling solutions for |e| < 1, a, (5 < 1, and a + /3 < 1 
. a)(l + ey- a (l - ef-? > 5. b) (1 + e) 1 ~ a (l - e) 1 "^ < 5, e > and (^b^) 2 """^! - 



a) 



b) 
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FIG. 7. Non trivial scaling solutions for |e| < 1 and a<l<(3ora<l<(3. 
a)(l + e) Q_1 (l - e)^ 1 > 1/5. b) (1 + e)^ 1 ^ - e)^ 1 < 1/5, and a < 1< /3 or a = 1, (3 > 1 and 
2^ > 1/5, or a < 1, /? = 1 and 2"- 1 < 1/(5. 




FIG. 8. Non trivial scaling solutions for |e| < 1 and a = (3 = 1 and 5 > 1. 
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